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Abstract 

Several VLBI EOP series were investigated with goal of determination of parameters of the Free Core 
Nutation (FCN). Both the amplitude and period of the FCN were studied using spectral and wavelet 
analysis. Our analysis reveals a variability of both the amplitude (known also from other investigations) 
and the period (or phase) of the FCN nutation. The FCN amplitude varies in the range about 0.1-0.3 
mas, and the FCN period — in the range about 415-490 solar days. The latter may be also explained by 
changes in the FCN phase. Comparison of time variations of the FCN period and amplitude obtained by 
different authors and methods shows substantial discrepancies at the edges of the period of observations. 

1 Introduction 

In this paper we investigate variability of the FCN parameters. Whereas variations of the FCN amplitude 
was already investigated (see e.g. [DH]), variations of the FCN period is not been studied yet. 

Modern theory of nutation predicts the steady FCN period of 431.2 sidereal days [3]. The FCN period 
also have been estimated from VLBI observations, and found to be about 430-431 sidereal days or about 
429-430 solar days (see, e.g. Table 4 in g]). 

In this paper we analyze four VLBI nutation series available in the IVS data base, sufficiently long 
and dense to obtain reliable estimates. We consider the differences between observed values of nutation 
angles and IAU2000A model (which is equivalent to MIIB2000 model without FCN contribution). For our 
purpose, we interpret the unpredicted part of observed nutation series in the FCN frequency band as the 
FCN contribution. 

2 Data used in analysis 

Four celestial pole offset series used in our analysis are BKG00003, GSF2002C, IAAO0201, USN2002B. We 
analyzed both raw (i.e. given on original epochs) and smoothed (equally spaced by 0.05 year) differences 
between observed nutation angles and the IAU2000A model. For smoothed series we also computed the 
weighted mean one. The parameter of smoothing was chosen in such a way to suppress oscillations with 
periods less then 1 month. Common time span for all series is 1984.0-2002.8. Figure [T] shows smoothed 
series used in our analysis. 

3 Analysis and results 
3.1 Spectral analysis 

For estimation of the power spectral density from both raw (unequally spaced) and smoothed (equally spaced) 
nutation series we used the Ferraz-Mello's method |5j which allows us to process both types of data. For 
supplement testing, we also compute the power spectral density using the Burg's method [6| . Figures [21— [4] 
show the normalized results of spectral estimation, and Table [T] presents the estimates of the FCN period. 
Comparison of results shows reasonable good agreement between the VLBI series, taking into account that 
we investigate rather week signal. 

The average estimated value of the FCN period is of about 434 solar days (about 435 sidereal days). 
This value is substantially greater than one found in j4j (431.0±0.6 sidereal days). However, when we used 
for spectral analysis only nutation series cut at the epoch 2000.2 which corresponds to the data span used 
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Figure 1: FCN contribution in the individual and mean series. 
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Figure 2: Spectra of raw data, Ferraz-Mello's method, solar days. 
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Figure 3: Spectra of smoothed data, Ferraz-Mello's method, solar days. 
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Figure 4: Spectra of smoothed data, Burg's method, solar days. 
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Table 1: Periods of the FCN contribution, solar days. 
Series Method BKG GSF lAA USN Mean 
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in [4], we obtain the FCN period of about 432 sidereal days which is close to found in [4] (see the last line 
in each section of Table [1]) . 

It is interesting to note a second period in the investigated frequency band, about 410 days. This period 
can be clearly seen in the Burg spectra, and as well it is discernible in the Ferraz-Mello spectra. Evidently 
a term with this period deserves further investigation^ 

3.2 Wavelet analysis 

At the next step we had applied the wavelet analysis to all the nutation series to investigate the time 
variations of the FCN period (phase) and amplitude, which is the main goal of our studjU. For this analysis 
we used program WWZ, developed by the American Association of Variable Star Observers and available 
as executable at the |http:/ /www. aavso.org/cdata/wwz.shtml Theoretical background of this method can 



be found in [7J . The results of the wavelet analysis are presented in Figures [SHU 

Figure [H] presents the final results of the present investigation. It should be mentioned that based on the 
comparison of FCN amplitudes found here and previous investigations [8] , we consider the results obtained 
before 1990 seems to be not very reliable. 

Of course, an important question arising from the obtained result is whether the variations of the period 
found from our analysis is an actual geophysical signal or an artifact caused by inadequate computational 
procedures. One can see that large increasing of the FCN period after «1998 corresponds to relatively low 
amplitude of the FCN oscillation. We have performed some tests to estimate how result of wavelet analysis 
depends on variable amplitude of input signal. 

For test purposes we used several artificial signals, and also we constructed new series as original one 
normalized by found variations of the FCN amplitude. The latter provides the FCN contribution series with 
near-unity amplitude with all other peculiarities inherited from real nutation series. The variations of the 
period found for the test series are practically the same as for real data. 

After all test, our conclusion is that found variations of the FCN period cannot be explained by compu- 
tational errors. Besides, the results of spectral analysis made for different subset of data also corroborate 
our conclusion. 



^After this paper was published in its original form, Krasinsky& Vasilyev (CMDA, 2006, v. 96, 219) gave a theoretical 
explanation of this FCN mode. 

^After this paper was published in its original form, the authors were pointed out that perhaps the first attempt to apply 
the wavelet analysis to investigation of the FCN was made by Schmidt & Schuh (ZfV, 1999, No. 1, 24). 
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Figure 5: Variations of the FCN period with time, raw data, solar days. 
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Figure 6: Variations of the FCN period with time, smoothed data, solar days. 
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Figure 7: Variations of the FCN amplitude with time, raw data. 
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Figure 8: Variations of the FCN amplitude with time, smoothed data. 
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Figure 9: Variations of the FCN period and amplitude with time; Atp * sin(e) (solid line), As (dashed line), 
and mean of and Ae (bold line). 
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One of the important points to be investigated is the edge effect which may lead to misinterpreting of 
the results of the wavelet analysis for the first and the last epochs. For this purpose we performed a special 
test with WWZ. We computed the FCN period and amplitude variations for several series starting with 
original one and cutting the first and the last 200 points from it. Then this process was repeated and in that 
way four test series were obtained, each starts 200 days later and ends 200 days earlier than previous one. 
Figures [TOl and [TT] show the result of this test. 

From this test, we can conclude that the edge effect may affect the result at the first and the last 500-1000 
days of the interval under investigation. Evidently, this fact should be accounted for to determinate the start 
point for prediction of the FCN contribution to nutation. 

4 Discussion and conclusions 

Variations of the FCN amplitude found in this investigation are close to ones used in the MHB2000 model, 
except the edge intervals (see Figure [T2|) . A reasons of these discrepancies may be insufficient quality of 
the VLBI data in the earlier 1980th and the edge effect present in the wavelet analysis results, as discussed 
above. 

As a supplement test, we computed variations of the FCN amplitudes immediately from the VLBI series. 
For this purpose we used the smoothed differences between observed nutation angles and model described 
above. 

Taking into account that the FCN contributions in Atp and Ae are two projections of the same variations 
in the Earth rotation velocity, we can compute the FCN amplitude as Amp{FCN) = y^d(A^^p+(i(Ae)^. 
Figure [13] shows the result of comparison. Comparing Figures [12] and [13] one can see that the direct com- 
putation of the FCN amplitude shows better agreement with the MHB2000 model. We can expect that 
amplitude estimates obtained with WWZ method are not good enough at the edges of the time interval (see 
also the discussion of the edge effect in wavelet analysis above). 

Comparison of the FCN phase variations found in this study and computed from the amplitudes of sine 
and cosine FCN terms of the MHB2000 model is presented in Figure [14] after removing the linear phase 
change corresponding to the FCN with permanent period. One can see that the FCN phase variations 
are similar in two approaches, though ours provides more smooth variations, and so for the FCN period 
variations. 

The FCN period (phase) most likely varies with time. Probably, change in the period is physically 
connected with change in amplitude. On the other hand, one can see that the variations of the FCN period 
show clear periodicity with a period about 5 years, whereas variations of the FCN amplitude does not show 
such an effect. 

Another reason of the observed behavior of the FCN period maybe variability of the FCN phase. Analo- 
gous effect was found also at the Chandler frequency [S], for which dependence of the period on amplitude, 
and the phase jump occurred during the period of the lowest amplitude were also found. 

It is interesting, that the Chandler wobble period also decreased in «1986-1988, and increased in wl989- 
1996 (see [ini[II])- Unfortunately, Polar Motion series studied in those papers are much shorter than one 
analyzed here to perform a reliable comparison. 

Variations of FCN amplitudes show several possible epochs of the excitation of the FCN, most of them 
are close to ones detected in [5]. 

Some tests we performed allow us to make a conclusion that investigated nutation series really contain 
such a signal with variable amplitude and period (phase). However, as stated above, we interpret the 
differences between observed nutation and the IAU2000A model as the FCN contribution, which may be 
too strong assumption. Possible interference of the FCN and other nutation frequencies should be carefully 
investigated. In particular, the authors of [T] pointed out a possible interference with near-yearly nutation 
terms, but the investigated period of observations seems to be long enough to separate these frequencies by 
spectral and wavelet analysis. 
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Figure 10: The edge effect in the wavelet analysis, the FCN period. 
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Figure 11: The edge effect in the wavelet analysis, the FCN amplitude. 



14 



FCN amplitude, mas 



0.4 



0.3 



0.2 



0.1 - 







/ s 



1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 

Figure 12: Variations of the FCN amplitude with time found in the present study (solid line) and implemented 
in the MHB2000 model (dashed line). 
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Figure 13: Variations of the FCN amplitude with time: solid line — original amplitudes, bold line 
smoothed amplitudes, bold dashed line — amplitude variations found from the wavelet analysis. 
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Figure 14: Variations of the FCN phase with time found in the present study (sohd line) and computed from 
the MHB2000 model (dashed hne). 
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